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Based on the time-dependent density-functional theory, we have derived a rigorous formula for 
the stopping power of an interacting electron gas for ions in the limit of low projectile velocities. 
If dynamical correlation between electrons is not taken into account, this formula recovers the 
corresponding stopping power of noninteracting electrons in an effective Kohn-Sham potential. The 
correlation effect, specifically the excitonic one in electron-hole pair excitations, however, is found 
to considerably enhance the stopping power for intermediately charged ions, bringing our theory 
into good agreement with experiment. 
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The concept of electron correlation permeates many 
subfields of condensed matter physics, but it is rather 
difficult to find good examples of its clear manifestation. 
It is even more difficult to unambiguously display its im- 
portance in time-dependent processes. Here we report 
a success in obtaining a clean theoretical formula di- 
rectly relating dynamical correlation, as characterized by 
the exchange-correlation (xc) kernel of time-dependent 
density- functional theory (TDDFT) , with the experimen- 
tally observed stopping power (SP) of solids for slow ions. 

The nonperturbative theory of scattering of noninter- 
acting electrons by a central potential combined with the 
density-functional theory (DFT) of Hohcnbcrg, Kohn, 
and Sham 1 ^ has proved successful in calculating the SP 
of an electron gas (EG) for slow ions, 3 This scheme accu- 
rately reproduces the measured energy loss of slow pro- 
tons and antiprotons in a variety of solids and success- 
fully accounts for the observed oscillatory behavior in the 
energy loss of slow ions with increasing charge, which is 
known to originate from the formation of closed shells 
of bound states around the ion^ These calculations were 
extended to projectile velocities approaching the Fermi 
velocity of the target but few attempts have been made 
to include the effect of many-body dynamical correlation 
that is absent in the current DFT schemes^ 

In this Communication, we address the effect of dy- 
namical correlation and we employ TDDFTZ. to derive a 
rigorous formula for the friction coefficient (SP divided by 
the projectile velocity at its zero value) of an interacting 
EG. This formula fully includes the effect of dynamical 
correlation through the imaginary part of a frequency- 
dependent xc kernel, and it reduces in the absence of 
dynamical correlation to the SP of noninteracting Kohn- 
Sham (KS) electrons first reported by Echenique et al£ 

By adoption of the widely used local-density approx- 
imation (LDA), our formula allows us to make a rea- 
sonably accurate evaluation of the friction coefficient of 



an interacting EG as a function of the projectile charge 
Z\e for Z\ less than about 15. Our calculations indicate 
that dynamical correlation considerably enhances the SP 
of slow ions, which is due to the excitonic effect in the 
electron-hole (e-h) pair excitations. In particular, this 
effect makes the ratio of maximum to minimum values in 
the Z\ -oscillatory behavior of the SP to decrease, leading 
to better agreement with experiment. 

We consider a recoiless probe particle of charge Z\e 
moving with velocity v in an otherwise homogeneous gas 
of interacting electrons at zero temperature. The average 
energy lost per unit length traveled by the probe particle 
(the so-called stopping power of the target) is simply the 
retarding force that the polarization charge distribution 
in the vicinity of the projectile exerts on the projectile 
itself. Accordingly, one can write (we use atomic units) 8 
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being the electron density induced by the projectile. 
Our starting point is the following rigorous expression 
for the friction coefficient of an interacting electron gas: 
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-AirZi/q 2 is the Fourier transform of the 
bare Coulomb interaction and Xq,k{w) is the double 
Fourier transform of the linear density-response function 
of a system of interacting electrons in the static field of an 
impurity of charge Z\ at r = 0. In order to prove Eq. J2J, 
we first expand the right-hand side of Eq. Q in pow- 
ers of Z\ and we then exploit the symmetry properties 
of many-fold density-response functions to demonstrate 
that this expansion divided by the projectile velocity at 
its zero value exactly coincides (to all orders in Z\) with 
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the corresponding expansion of the friction coefficient of 
Eq. 

The density-response function Xq,k(w) depends on the 
unknown ground and excited states of a many-electron 
system and is, therefore, difficult to calculate. Neverthe- 
less, TDDFT allows us to express Xq,k( w ) as a solution 
of the integral equation^ 



Xq,kM= / rfpe q ,p(w)Xp,k(^ 
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where e qk (^) represents the so-called inverse test- 
charge-electron dielectric function, 

q]»^ q , k -/ dpx q ,p(^)[(4^ 2 )^,k+/^k(^)] , 
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and / qk (w) is the Fourier transform of the xc kernel 
f >o](r,r» = ^>](r,c,)/<Kr>)U„ , (5) 

with v xc {n](r,uj) being the frequency-dependent xc po- 
tential. Xqk^) ^ s ^ ne Fourier transform of the density- 
response function of noninteracting KS electrons, i.e., in- 
dependent electrons in the effective KS potential 

v(r) = -Zx/r + [ dr'n (r')/|r - r'| + v xc [n ](r), (6) 



with v xc [no] (r) being the static xc potential at the elec- 
tron density no(r). The unperturbed density no(r) is 
obtained by solving self-consistently the KS equations 
of DFT with the effective potential v(r) of Eq. JBJ. At 
u) = 0, we hnd the following sum rules: 
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where n q and Vq are the Fourier transforms of the density 
n (r) and the potential v(r). With the aid of these sum 
rules, substitution of Eq. © into Eq. yields 



Q = Qi + Q 2 , 
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Equations l|%)l-p()[l. which describe the SP of an in- 
teracting many-electron system in the low- velocity limit, 
show that the low-velocity SP of interacting electrons 
[Eq. |[5J|] can be obtained as the sum of (i) the SP of 
noninteracting KS electrons [Eq. (JSJ)] and (ii) an excitonic 
contribution [Eq. QIC)))], which accounts for the dynamical 
correlation that is absent in Eq. (JSJ. 



In order to show that the SP of Eq. @ is precisely the 
SP reported in Ref . under the assumption of individual 
elastic electron scattering by the effective KS potential of 
Eq. © , we write an explicit expression for the imaginary 
part of X q k( w ) ty using incoming and outgoing single- 
particle states, Ip'^) and |p^ + ^), and the corresponding 
single-particle energies, e p < and Ej/Q 

Im Xq>k M = (2tt)- 2 J dpdp'l.f(e p ,)-f(e p )}6(Lo-e p +e p ,) 

x (p' ( ~V qr |P (+) ><P (+) |e~ ikr |p' (_) >> (11) 

where f(e) denotes the Fermi distribution function. Sub- 
stituting Eq. into Eq. (JHJl and referring to the Lipp- 
man-Schwinger equation of scattering theory, one finds 
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where k F is the Fermi momentum, o- tr (q) is the so-called 
transport cross section, and n is the density of the ho- 
mogeneous EG in the absence of the projectile. 

An evaluation of Eq. © [or, equivalently, Eq. (T2J] 
and Eq. IjlOl) involves the unknown quantities i> xc [rio](r) 
and f xc [no](r, r'; u>), respectively. Existing DFT calcula- 
tions from Eq. I|12|) have been implemented by invoking 
the LDA, i.e., by replacing the xc potential v xc [no](r) 
by the xc potential of a homogeneous EG with the lo- 
cal density no(r). At the same level of approximation, 
we use the LDA to evaluate the excitonic contribution 
of Eq. (|10l) . In this approach, the frequency-dependent 
(nonadiabatic) xc kernel f xc {no](r, r'; ui) is obtained asi£ 
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where f^ c (n,u>) denotes the q — > limit of the frequency- 
dependent xc kernel of a homogeneous EG of density n. 
Due to the spherical symmetry of Uo(r), substitution of 
Eq. O into Eq. (TDJ yields 
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An explicit expression for the imaginary part of 
f£ c (n,u) has been reported recently, which is exact 
at low frequencies to leading order in the Coulomb 
interaction*^ Since the frequency derivative entering 
Eq. (|14fl is negative definite,-- this equation shows that 
dynamical correlation enhances the SP of a system of in- 
teracting electrons. This enhancement can be ascribed to 
the attractive e-h multiple scattering or excitonic effect 
in the e-h pair excitation^ 

The calculations presented below have been carried 
out from Eqs. lfT2|l and |TI)l. with use of the LDA xc 
potential v xc [no](r) of Perdew and ZungeoSil and the 
frequency-dependent LDA xc kernel f%°(n,uj) of Qian 
and Vignaleii^ We have performed calculations for var- 
ious values of Z\ and the electron-density parameter 
r s = (3/47rn) 1 / 3 . We have found that dynamical cor- 
relation can be neglected at small projectile charges and 
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FIG. 1: Friction coefficient of a homogeneous EG versus the projectiie charge Z\ (a)-(c) and the EG-density parameter r s (d), as 
obtained in the LDA from Eq. 1121 (dashed chained curves) and from the sum of E as. (I12H and 1141 (solid chained curves), (a) 
Solid circles, triangles, and diamonds are transmission measurements from Refs.llll.ll21 and[l3jof the random SP of C for ions 
moving with velocities v = 0.41, 0.83, and 0.25a.u., respectively, (b) Solid circles and triangles are transmission measurements 
from Refs. |n] and ^3 of the random SP of Al for projectiles moving with velocities v = 0.41 and 0.83a.u., respectively. Solid 
squares (stars in the inset) are the measured SP of Al of Ref. Ilia) for slow ions (v = 0.5a.u.) moving at a distance of 1.0 
(1.2) a.u. from the last atomic row at the surface, (c) Solid squares are the measurements of Ref. of the SP of Au for slow 
ions (y — 0.68a.u.) channeled along the (110) direction. These measurements are compared to EG calculations with r s = 2. 
The dotted line is the calculation of Ref. |6J with dynamical xc effects included in the framework of linear-response theory. 



high electron densities, where many-body effects play a 
minor role. At metallic densities (1.5 < r s < 6), however, 
the correlation effect is found to considerably enhance the 
SP for intermediately charged ions. 

At low velocities, the energy loss of ions in solid ma- 
terials is mainly due to the stopping power of valence 
electrons. Since these electrons are known to be well de- 
scribed by a homogeneous EGj2i we apply our theory to 
the cases of C, Al, and Au, and compare the result of 
our EG calculations to the available experimental data 
for these materials (see Fig. 1). 

First, we compare the results we have obtained for the 
friction coefficient of a homogeneous EG with r s = 1.59 
and r s — 2.07 [solid chained curves in Figs. 1(a) and 
1(b)] to the transmission measurements of the random 
SP for slow ions in C and Al (solid circles, squares and 
diamonds). If we neglect dynamical correlation (dashed 
chained curves), our calculations reproduce those of 
Ref. 4. The effect of dynamical correlation is found to be 
negligible at low ions charges; however, as Z\ increases 
this effect becomes comparable to the SP of noninteract- 
ing KS electrons. We observe that the impact of dynami- 
cal correlation becomes relatively larger at the formation 
of closed shells of bound states around the probe particle, 
where the SP has a local minimum. As a result, the ratio 
of the maximum to minimum values in the oscillations of 
the SP decreases in the presence of dynamical correlation, 
leading to better agreement with experiment. 

Differences between our calculation and the transmis- 
sion measurements of the random SP of C and Al should 
be partially due to the contribution from elastic colli- 
sions of the projectile with the atomic cores of the tar- 
get, which is not included in our calculation. Hence, we 
have also considered the measured SP of Al for slow ions 
moving at a distance of d =1.0 a.u. from the last atomic 



row at the surface^ [solid squares in Fig. 1(b)], which 
is free from elastic collisions with target atoms. At this 
distance, the electron density is close to the bulk value 
(r s = 2.07) and the SP is expected to be only slightly 
smaller than that for ions moving in a uniform EG. 22 Fig- 
ure 1(b) shows that this is indeed the case, although our 
many-body calculation (solid chained curve) is too high 
for Z\ > 15. As Z\ increases the electron-density vari- 
ation near the projectile becomes large, and differences 
between our calculations and the experimental data may, 
therefore, be ascribed to the failure of the LDA. 

To elucidate the effect of the inhomogeneity of the EG 
near the Al surface, in the inset of Fig. 1(b) we plot 
again our many-body calculation at r s — 2.07 and the 
measured data at d =1.0 a.u., but now together with the 
measured SP of Al for slow ions moving at d =1.2 a.u>i£ 
(solid stars) and the calculated curve for the correspond- 
ing density of r s — 2.2 (solid chained curve with stars). 
Since the two theoretical curves essentially coincide, this 
figure shows that at 1.2 a.u., where the gradient of the 
electron density is considerably larger than that at d =1.0 
a.u.ji&i SP is worsly described by the homogeneous EG. 

Contributions from elastic collisions with the target 
atoms should also be absent in channeling experiments, 
in which the projectile penetrates a solid along a high- 
symmetry direction. In Fig. 1(c) we compare the ex- 
perimental data of Ref. (solid squares) for slow ions 
channeled along the (110) direction in Au with the result 
of our many-body calculation for the friction coefficient 
of a homogeneous EG with r s = 2.0. This figure shows 
that many-body correlation effects tend to decrease the 
ratio of the maximum to minimum values in the oscilla- 
tory behavior of the SP, thus bringing the SP into nice 
agreement with experiment for projectile charges up to 
Z\ =19. However, comparison between theory and ex- 
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periment at Z\ > 19 suggests that at these large values 
of Zi the strong inhomogeneity of the EG suppresses dy- 
namical correlations and calls for going beyond the LDA. 

Finally, Fig. 1(d) displays the interacting and non- 
interacting SP for H + and He ++ , as a function of r s . 
Figure 1(d) shows that (i) the effect of dynamical corre- 
lation is very small in the case of H+ , in which case the 
electron density variation is comparatively small, and (ii) 
the dynamical correlation effect becomes comparatively 
more important as the electron density decreases. 

In conclusion, we have presented a rigorous TDDFT 
for the stopping power of an interacting EG in the limit 
of low projectile velocities. In the absence of dynami- 
cal correlation, our theory yields the SP of noninteract- 
ing KS electrons, which we have shown to be equivalent 



to the existing DFT-based elastic electron scattering ap- 
proach. We have investigated the effect of many-body 
dynamical correlation by using an accurate representa- 
tion of the long-wavelength frequency-dependent xc ker- 
nel. Our results show that dynamical correlations con- 
siderably enhance the SP of intermediately charged ions, 
thereby bringing the theory into better agreement with 
the experimentally observed oscillations of the SP for pro- 
jectiles of charge Z\ <15. However, at higher values of 
Z-y, one needs to go beyond the LDA in the description 
of the dynamical xc effects. 
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